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C^ ■ Abstract 

In this paper we give a precise estimate of the discrepancy of a 
CN ' class of uniformly distributed sequences of partitions. Among them 

we found a large class having low discrepancy (which means of order 
j^). One of them is the Kakutani-Fibonacci sequence. 
< 

u 

^ : 1 Introduction and preliminaries 

-)— » 
a . 

S. Kakutani studied in [K] the notion of uniformly distributed sequences of 
partitions of the interval [0, 1]. He introduced the following construction. Fix 
a number a G ]0, 1[. If vr is any partition of [0, 1], its a-refinement, denoted 
>• ' by avr, is obtained subdividing the longest interval(s) of vr in proportion 

^ ■ Ci/i^ ~ <^)- By a^Ti we denote the a-refinement of a"~^7r. 

O ■ Let a; = {[0, 1]} be the trivial partition of [0, 1]. The sequence {a'^u} will 

^ ! be called the Kakutani a-sequence. 

in 

o 
o 



Definition 1.1 Given a sequence of partitions {vr„} of [0, 1], with 



nn = {[yUyn^<^<Hn)}, 

r> I we say that it is uniformly distributed (u.d.) if for any continuous function 

c^ I /on [0, 1] we have 



k(n) 



"^°o k{n) 



i=l 



]^:SLTr:r^T.f(^y'l)= L fit)dt. (i) 



We can now state Kakutani's result. 
Theorem 1.2 The sequence {a"'U}} is uniformly distributed. 



This result caught the attention of several authors in the late seventies, 
when a different proof of the theorem was given ([AF]), and other papers 
were devoted to a stochastic version of it (see [LI], [L2], [PvZ], [vZ]). 

Recently, the procedure introduced by Kakutani has been generalized in 
several directions. 

In [ChV] this notion has been extended to compact metric spaces. 

In [V2] Kakutani's splitting procedure has been generalized, producing a 
new class of u.d. sequences of partitions as follows. 

Definition 1.3 Consider any non trivial finite partition p of [0, 1]. The p- 
refinement of a partition tt of [0, 1] (which will be denoted by pvr) is obtained 
by subdividing all the intervals of it having maximal length positively (or 
directly) homothetically to p. 

Obviously, if p = {[0, a[, [a, 1]}, then the p-refinement is just Kakutani's 
a- refinement. 

As in Kakutani's case, we can iterate the splitting procedure. The p- 
refinement of pir will be denoted by p^vr and, for all n G IN, the p-refinement 
of p"~^7r will be denoted by p""K. If {p"'uj} denotes the sequence of successive 
p-refinements of the trivial partition u, the following theorem holds (see [V2, 
Theorem 2.7]): 

Theorem 1.4 The sequence {p"cu} is uniformly distributed. 

In [CVl] the splitting procedure has been generalized to higher dimen- 
sions, providing a sequence of nodes in the hypercube [0, l]'^ which is uni- 
formly distributed. 

In [IV] the authors study uniform distribution on fractals. 

In [CV2] the authors present a von Neumann type theorem (see [vN]), 
which provides uniformly distributed sequences of partitions of the interval 
[0, 1] starting from sequences of partitions whose diameter tends to when 
n — 7- oo. 

In [DI] the authors give upper estimates of the discrepancy of p-refinements 
of the interval [0, 1]. 

Actually, the theory of uniformly distributed sequences of partition is 
deeply connected to the theory of uniformly distributed of sequences of 
points, which have a large application in quasi-Monte Carlo methods (see 
[N]). 



For a large overview on uniformly distributed sequences of points, see 
[DT] and [KN]. 

Let us consider now a special class of p-refinements of the trivial partition 

u. 

Definition 1.5 Let us fix two positive integers L and S and let < /3 < 1 
be the real number such that L/9 + S/3^ = 1. Denote by p l^s the partition 
defined by L "long" intervals having length /3 followed by S "short" intervals 
having length /3^. By {p^s^} i'^^ (Pls) ^^ short) we denote the sequence 
of successive pL^s-refinements of the trivial partition u. They will be called 
LS-sequences. 

It is clear that the partition p^ s is obtained by dividing all the longest 
intervals of p^l~g homothetically with respect to Pl,Sj and that each partition 
p2 g contains only two kinds of intervals: the long intervals have length /?"■ 
while the short ones have length /S""*"^. 

By Theorem 1.4 we know that the sequence {p^ 5} is uniformly dis- 
tributed, but for this particular simple class of sequences a direct proof was 
given in the lecture notes of a course on Uniform Distribution delivered by 
A. Volcic during the "Workshop on Measure Theory and Real Analysis" held 
in 2002 ([VI]). 

In this paper we give precise estimates of the discrepancy of LS*- sequences. 

Among the family of LS'-sequences having low discrepancy, there is the 
Kakutani a-refinement with a + a"^ = 1 (or, equivalently, L = S = 1). 
Apart from the obvious case a = 1/2, this is the only Kalutani sequence for 
which the exact discrepancy is now known. This sequence may be called the 
Kakutani-Fibonacci sequence. 

2 L/S-sequences of partitions 

This section will be devoted to the study of LS'-sequences of partitions and 
their discrepancy. 

Denote with tn the total number of intervals of p^ 5, with /„ the number 
of its long intervals and with s„ the number of its short intervals. Obviously 
we have the following relations: 

^n ^n ~T~ f^n7 'n -^ ^n—1 ~r ^n— 1; ^n *^ 'n— 1* 



The terms of the sequence t„ can be calculated solving the following 
difference equation 

tn ^ -Lj tn~l ~\~ »-' tn-2 

with initial condition to = 1 and ti = L + S . 
We easily get 

We note that /„ and s„ satisfy the same difference equation, but the initial 
conditions are different since Iq = 1, li = L and Sq = 0, Si = S, respectively. 
It is convenient to introduce the constants 

, 1 + 5/3 , ^ S/3 - S/32 

which will frequently appear in the sequel, so that formula (2) can be rewrit- 
ten in the following way 

_ A- B{-S/3'r .o^ 

~ P^ ■ ^^ 

Since 5/3^ < 1, we see that t„ has the same order as -^ when n — )■ 00. 

The following remark will play a role in identifying the different cases in 
the statement of Theorem 2.3. 

Remark 2.1 From L(3 + 5/3^ = 1 we get 



_ VL2 + 45 - L 
' 25 

from which we derive that 

5/3 < 1 if and only if 5 < L + 1, 
5/3 = 1 if and only if 5 = L + 1, 
5/3 > 1 if and only if 5 > L + 1. 

Now we recall the definition of discrepancy and star- discrepancy. 



Definition 2.2 Given a finite subset W = {wi,W2, ■ ■ ■ , ^Af} of tlie interval 
[0, 1], tlie discrepancy of W is defined as 



D{W) = sup 

0<a<fe<l 



N 



j;rJ2X[a,b[iuJj)-{b-a) 



while the star- discrepancy is defined as 



D^'iW) = sup 

0<fe<l 



N 



N 



J2x[o,b[{wj)-b 



i=i 



If we consider a sequence of partitions {7r„} of the interval [0, 1], with 
^n = {[y"_i,|/r[, 1 < "^ < k{n)}, we may define 



D{ 



TT^ 



sup 

0<a<b<l 



k{n) 



k{n) . 1 



Exkm(^7)-(& 



and 



^*(7r„) = sup 

0<6<1 



k{n) 



k{n) fr{ 



Ex[o,m(^")-^ 



(4) 



It is well known that {vr^} is uniformly distributed is and only if D{'Kn) -^ 
when n ^ oo, and that D*{W) < D{W) < 2D*{W). 

The speed of the discrepancy is important in the applications. It is also 
well known that the best discrepancy is - and it is attained, for example, by 
Knapowski's sequence {[^--^, -[,!<?< n\. 

We are now ready to determine the star-discrepancy of the sequence 

{Pl,s}- 

Theorem 2.3 The discrepancy of the sequence of partitions {p^s)' O'S n ^ 
oo, behaves in the following way: 






m) D{pls) -J^^fS>L + l, where 7 = 1 + '-^. 



Proof. We denote the intervals of p^ s ^^ follows 

Pls = M-i,ynA<^<tn)} 

and we study the behavior of the star-discrepancy of p2 5- 

Denote by Lp and Sp, respectively, a long interval and a short interval of 
the partition p^ g, and first of all let us evaluate for n > p the differences 

and 

rT.xsM)->^(Sp), 

where of course A denotes the Lebesgue measure and A(Lp) = (3^ and \{Sp) = 
(3P+\ 

We observe that, since pj g is a refinement of p^ 5, each long interval Lp 
of p^ g is the union of consecutive intervals of pj 5 (with n > p) and, since 
the splitting procedure is "self-similar" , this union reproduces up to a factor 
/3~^ the partition p2^- Therefore we have, for any n > p, 

J2XLpiyJ) =tn-p- 

i=i 
With a simple calculation we obtain, for n > p, 

= -- {-spy [{-spY'p - pp] . (5) 

^n 

If we consider now a short interval Sp of p^, 5, it becames a long one in 
P^L^Si which means that it is of the type Lp+i, so we can rewrite formula (5) 
in the following way: 



in j=i 



^X5,(y;)-A(5,) = ^^if^-/3^ 



= -- (-5/3)" [{-S^)-P-^ - Z^^'+^J (6) 

for every n > p + 1. 

Fix now b e]0, 1[. 

Let [6"~^,62~^[ be the interval of p^'^g containing b. Since the number of 
points of p^~g contained in [&i~^, &2^^[ is 1 if this interval is a short one, and 
it is L + S* if it is a long one, we have 

1 in 1 in 

-rT.xioMiy^)-b<-j:x[o,,r^M)-br' = 

1 tn r -I- S* 



tn j=i 



and 



1 in 1 in 






1 tjT, 1 Cn 



1 in r -I- t? 

= rEx[o,.-[(^;)-^r'-^- (8) 

We will estimate (4) using the upper and the lower bounds in (7) and 
(8), so that we will evaluate j- Z]j=i X[o h^-'^i iVj') ~ K~^ f°^ k = 1,2. To this 
purpose it will be convenient to represent [0, b^~^[ (from now on we will write 
[0, 6"~^[ for short) as the union of consecutive intervals stemming from the 
partitions p^g ior i < n — 1. 

Consider first all the consecutive intervals ll,!^, . . . , I^^ of p\^ g such that 
{JT=i 1} C [0, 6"~H (mi could be zero). 

Next, take all the consecutive intervals If , /|, . . . , /^^ of Pl s contained 
in [0, 6"~^[ \ lJi=^i 1} (again m2 could be zero). 



Proceed this way taking, at the final step, all the rrin-i (possibly zero) 
consecutive intervals /"~^, l2^^, . . . , -^m^^i ^^ Pls ^^'^'^ ^^^^ 

'm„-i n-2 /rrip \ 

U ir' c [o,6-M\ u U^fl 

i=l p=l \i=l / 

SO that at the end 

[0,b-'[= U U^f • (9) 

Thus [0, 6"^^[ is represented by the union of (at most) n — 1 blocks of 
consecutive intervals of p^ g for 1 < p < n — 1. 

Let /p and s^ be, respectively, the number of long intervals (denoted by 
Lf) and short intervals (denoted by Sf), respectively, of p^, ^ contained in 
[0, b[ (obviously, Ip + Sp = nip). Then using (5), (6) and (9), in order to 
evaluate the star-discrepancy we can write: 

-^ tn /n-l /mp ^ 

= ^g'^(u;;;(u»'.'.-))<^'>-MUiy/f^ 



j=i y^p^i\^i=i ijj \p=i \i=i 

n-l nip / -r t„ 
p=l j=l y-n j=i 
n-l ^p I T t„ \ n-l ^v I -[ tn 

E E - E XLfd/;) - A(LD + E E - E x^r (y,") - A(5f ) 
p=i i=i y'-n j=i y p=i 1=1 y*-" j=i 

R n-l 

= —(-5/3)" E {-? [(-^/3)"' - /?'] - 4 [(-^/3)"'"' - /3"+']} ^ 

^n p=i 

D n— 1 

= —{-SPY E f-^J [(-^/3)"" - /31 - 4 [(-^/3)"""' - /?"+'] 



f-n 



p=l,p even 



R n— 1 

+—(-5/3)" E Hp[(-^/3)~'-/31-sM(-5/3)-^-^-/3^+^] 



" p=l,p odd 



B 



n— 1 



n-1 



^p=l,pevcn p=l,p odd 



rj I n— 1 n— 1 

^ I p=l,p even P=1,P odd 



(10) 



It is clear from the construction that ntp <L + S — l,0<lp<L and 
< Sp < aS* for all 1 < p < n, otherwise at least one of the intervals of lj™f^ If 

would be already present in USi^ ^f^ ■ 

We have to distinguish the cases when S/3 ^ 1 and SP> = 1. 

We will give upper and lower estimates of (10) in the case n even and in 
the case n odd. 

Assume first that n is even and let n = 2h, h > 1. 

We note that the first term in (10) is positive, and the second term in 
(10) is negative since B > and {Sl3)~'^ + /3^ = "gia • Than with simple 
calculations we get the following upper and lower estimates: 



.B{L+S-l)^-^^l2 



^2h 



1 - {SP) 



2h 



2h 



1-/3 

(l-52^2)(^^)2/.-2 1_^2 



+ 



(S^) 



2h 



p-zny < 



r- 



^2h j=i 



\2h 



t2h I 



1 - {SP) 



2h-2 q2 _ o2h 

;i - 52/32) (5^)2/.-! + ^ 1 _ /32 



5/3 



/3 



:ll) 



Assume now that n is odd and let n = 2h+l, with /i > 1. Then with the 
same arguments from (10) we derive the following upper and lower bounds: 



_S(i + S-l)M!^x 



hh^ 



X l2 



1 - (5/3) 



2h-2 



:i -52/32) (5/3)2^- 



-1 ^ 1-/32 



J_ , « I ^ I «2fe+l 

5/3 '^ (5/3)2^H 



T + /5^ 



< 





1 t2h+l 

hh+i j=i 


X[o, 


,..^{yr')-b'' 


< 




<2B{L + S 


hh+i 


] 


I - (5/3)2^+2 


1- 


■ /32^+2 


(1 


-s^p^){spy'' 


1 


-P' 


)S/3 = 1 












equation (10) 


becomes 












1^ 


,6"- 


4y-)-b-' = 







(12) 



p I n— 1 n— 1 

'^ (_p=l,p even P=1,P odd 

75 I n— 1 71—1 I 

-— (-ir E ^p,l[1-/31+ E m,,,[l-/3^+i] . (13) 

" ^_p=l,p even P=1,P odd j 

If we assume n = 2h, with h > 1, from equation (13) we derive the 
following upper and lower bounds 



BiL + S-1) i^^ , ^ 1-/3 



2h 



2/1 + 1-2 ^ - /3'^ < 



t2?. I 1-/3 

'-2/1 j=i 

If n = 2/i + 1, with h > 1, from (13) we derive the following inequalities 
t2h+l [ 1-/3^ J 

1 *2h+l 

<r^Ex[o,.-[(yr^)-&^^< 
10 



^B(L + S-l)f _2/3=-r« 



t 



2/1+1 



1-/32 



(15) 



Since we enclosed the generic 6 G [0, 1[ into the interval [&i^^,&2^^[ and 
we obtained the corresponding estimates (11), (12), (14) and (15), the same 
estimates hold for 6 G [0, 1[ and therefore they can be applied to evaluate 
the star discrepancy. 

It is time now to distinguish among the three cases explicited in the 
statement of our theorem, in order to prove the various asymptotic behaviors 
of the star-discrepancy. 

i) S <L + 1 

From Remark 2.1 this condition is equivalent to the case Sf3 < 1. 
If n is even, from (7), (8) and (11) simple calculations give 



X l2 





~[L + S + B{L + S-l)x 






1-52^2 y^P) l_/32^ 


+ 1 


- {sp^f^ 


< 


<t2h 


sup , } ,x%b[iyf) 

0<6<1 \J'2h j=i 


-b 
1 . 


< 





X <2 



1 - S^/3^ 



<L + S + B{L + S-l)x 



^^-i^^r-\^^,^y.Pj-P- 



P' 



(16) 



If n is odd, from (7), (8) and (12) we obtain 

^2/32 _ ^s(5f'' 



-[L + S + B{L + S-l)]l2 
-[L + S + B{L + S-l)] 

< ^2/1+1 



02 _ o2h 
2fe+l ^ ^ 

1-52/32 ' ^^^^^ 1-/32 



+ /3{S/3) 



2h 



+ /3(5/3)2'*+i + l + (5/32)2''+!} 



< 



1 '-2h + l 

0<6<1 \ t2h+l j^i 



sup 



2h+l\ 



-h 



< 



11 



<L + S + B{L + S-1){2 



Sf3 - (g/3)^^+3 
1 - 52/32 



(5/3) 



2/1+1 



1-/3 



2h+2 



1-/32 



■ (17) 



Since /3 < 1, S*/? < 1 and S/S'^ < 1, there exist the following limits: 



lim {2 



5/3 - (5/3)2^-1 , «,em2h/3'-/3'' 



1 - 52/32 



+ /3(5/3)^ 



1-/32 



+ (5/3)2'*-^ + /3(5/3) 



2h 



lim <^2 



5/3 - (5/3)2 ^+3 
1 - 52/32 



-(5/3) 



2?i- 



-1 1-/3 



2/1+2 



1-/32 



25/3 



1 - 52/32 



and 



lim {2 



52/32 - (5/3)2^^ , ..o.^2/.+i/3'-/3'' 



1-52/3 



2«2 



+ /3(5/3)^ 



1-/32 



+ 



+{SPr + /3(5/3)2'*+i + 1 + (5/32)2'*+!} 



lim <^2 



A— >oo 



52/32 - (5/3)2^^+2 

1 - 52/32 



(5/3) 



1 - /32'' 

1-/32 



+ 1 - (5/32) 



2\2/i 



252/3 



2/92 



1 - 52/3^ 



+ 1, 



SO that the sequences of upper and lower bounds in (16) and (17) are bounded. 
Then we conclude that there exists a positive constant Ci, independent on 
n, depending only on L and 5 (actually, /3 depends on L and 5), such that 
for all n e IV 



tn sup 
0<6<1 



1 tn 



n j=i 



<Ci, 



which implies that, as n — )■ 00, 



12 



D*{{pls}) - -. 



ii) S = L + 1 



First of all we note that, by Remark 2.1, we have Sf3 = 1. 
If n is even, from (7), (8) and (14) we obtain 



B{L + S-l] 



< 



log t2h 



1 - /92 



1-/3 



2h 



< 



logt 



2h 



1 fr'' 

sup I ^J2x[o,b[iy]'')-b 

0<b<l \ t2h j=i 



< 



<B(^+5-i)L,_,^,^g!^^; 



log t2h 



1-/3^ 



(18) 



If n is odd, we obtain from (7), (8) and (15) that 



B{L + S-1) [^, , ^^13^-/3^'' 



< 



log t2h+l 
t2h+l 



2/1 + 2/3- 



1-/32 



+ /3 + /^'^^^^ < 



logt 



2/1+1 



sup 7 — T.XioMiyT^')-b 



0<b<l \ t2h+l j=i 



< 



< B(L + S-l) l_^p-'-p-'» 



logt2h 



+1 



1-/32 



(19) 



From formula (3) we have 



logt„ = log(^^-i?(-5/3^ 



log 



'A-B{-S(3 



2\n' 



\og{A - B{-S/3T) - nhg /3 



13 



for all n G W so that, since /3 < 1 and SfS"^ < 1, 



lim 



2h-l + 2f3^ + l3 

log t2h 



1-/32 



lim 

h—¥oo 



2/1+1-2^-/3 

log t2h 



nh 



2/1-2 



o2_o2h+2 

r^;s2 



lim 

h^oo log t2/i+l 



lim 



2/1 + 2/3^3;:^ + /3 + /3 



2h 



logt 



2h+l 



log/9 



Taking the above limits into account, we conclude from (18) and (19) 
that there exists a positive constant C2, independent on n, depending only 
on L and S, such that 



1 — T *^^P 

login 0<b<l 



1 tn 



<C2, 



which implies that, as n ^ 00, 



D*{{pls}) 



logt^ 



iii) S> L + 1 

By Remark 2.1 this condition is equivalent to S/S > 1. 
If n is even, from (7), (8) and (11) we have: 



B{L + S-l] 



/3^H: 



2h 



1 - (5/3) 



2h 



2h 



'l-S^/3^){S/3) 



2/1-2 



1-/3 
1-/32 



+ 



(5/3) 



2h 



- /3''' < 



< 



(5/3 



2\2/i 



s^P :;— IIx[o,m(^j )-^ 



0<b<l \ t2h ,=1 



< 



^B(L + S-2)L 



/32'»t 



2/1 



1 - (5/3) 



2/1-2 



+ /3 



/32 - /32^ 



(1- 52/32) (5/3)2''-i "^ 1-/3 



^+4. (20) 



If n is odd, from (7), (8) and (12) we get 



14 



B{L + S-l] 



X 



X {2 



1 - {SP) 



2h-2 q2 _ Q2h 

[1 - 52/32) (5/3)2/^-1 + ^1-/32 



^^ + /^ + (^/3)2/.+l + ^ 



2fe+l 



< 



< 



1 



(sp 



2\2h+l 



1 t2h+l 

sup I Y. XloMiV] 

0<b<l \ 12/1+1 j=i 



2/.+U _ ^ 



< 



< 



2E(L + g-l) 



1 - (5/3) 



2/1+2 



1-/3 



2/1+2 



;i-52^2)(^^) 



2/i 



1-/32 



(21) 



Since S/3 > 1, Sf3'^ < 1 and /3 < 1, using formula (3) we have 



1 



/i^oo /32''t2/» 



1 - (5/3) 



:i- 52/32) (5/3) 



2/1-2 o2 _ o2h 

2/.-1 ^1' 1 _ ^2 



+ 5^ + ^ 



lim 



X 



X <^2 



1 _ (g/3)2fe~2 /32-/32fe 

:i- 52/32) (5/3)2^^"! +^ 1-/32 



+ ^ + ^. 



lim 

/,^oo /32'^+lt2/^+l 



{2 



1 - (S;3) 



2/1-2 



:i- 52/32) (5/3)2'* 
1 



^ ^^_^2^ 
'1 ~^^ 1 _ /^2 



-/32 



+ 



5/3 ^ (5/3)2''+i ^ '^ ^ 



^^ 



1 /33 



5/3(52/32-1) 1-/32 



5/3 



/3 



and 



15 




/3 



2h 



h^^ /32^+lt2h+ 



1 



From (20), (21) and the previous limits, we conclude that there exists a 
positive constant C3, independent on n, depending only on L and S, such 

that 



1 



[Sp^P 0<6<1 



1 tn 

-rJ2x[0MiyJ)-b 



<c. 



for all n G W, which implies that, as n — )> 00, 

^({p2,5}) ~ {sp'Y 

It remains only to observe how (S*/?^)" can be written in terms of t„. 
In fact, {S(5^Y = /3t, where 7 = 1 + ^^|^ < 1 since S(5^ = {S(5)(5 = (5^+^ 

with c = °^^ f -* < 0. Consequently, (S*/?^)" and ^7 have the same order at 
infinity since 



2\n 



l/{tnY 



ip 



"17 



[A + B{-Sf3^Y 






[A + B{-Sf3^Y 



^A^ 



as n — )■ 00, and finally, as n — )■ 00, we have 
where 7=1 



1 



log(g/3) 
log/3 • 

Since the discrepancy and the star-discrepancy are equivalent, the theo- 
rem is proved, -k 



16 



3 Conclusions 

Uniformly distributed sequences of partitions having low discrepancy turn 
out to be useful in quasi-Monte Carlo methods, as it is clear from formula 

(1). 

LS'-sequences of partitions offer the advantage that each partition is the 

refinement of the previous one. Nevertheless, the disadvantage of these se- 
quences lies in the fact that the function has to be evaluated only on "blocks" 
of points, so that it is not possible to use intermediate number of points be- 
tween tn and tn+i- 

There exists a connection between u.d. sequences of partitions and u.d. 
sequences of points. In [V2, Theorem 3.4] it has been proved that a random 
reordering of the points of a u.d. sequence of partitions gives with probability 
one a u.d. sequence of points. That theorem, even if interesting from a 
theoretical point of view, does not give any information about the discrepancy 
of these random sequences of points. 

The author conjectures that it is possible to associate to LS'-sequences 
having low discrepancy sequences of points having low discrepancy too, i.e. 
of order -^^jj—- Some encouraging results have been already obtained, and will 
be presented in a forthcoming paper. 
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